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We discuss an experimentally amenable class of two-particle states of motion giving rise to nonlocal 
spatial interference under position measurements. Using the concept of modular variables, we derive 
a separability criterion which is violated by these non-Gaussian states. While we focus on the free 
motion of material particles, the presented results are valid for any pair of canonically conjugate 
continuous variable observables and should apply to a variety of bipartite interference phenomena. 
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Is it possible to deduce entanglement from an inter- 
ference pattern, to perform an "entangled Young exper- 
iment," following the famous single-particle interference 
experiments? While the wave-particle duality of single 
material particles has been a central theme since the early 
days of quantum mechanics and is impressively confirmed 
in interference experiments [H @] i a similarly convincing 
demonstration of quantum nonlocality as implied by en- 
tanglement has proven to be much harder to implement 
with matter waves. 

Although recent experimental progress, in particular 
in controlling ultracold atoms, has rendered experiments 
conceivable that probe entanglement in the free motion 
of material particles, a direct implementation of most 
schemes that have proven successful with other contin- 
uous variable degrees of freedom (e.g., field modes) fails 
due to the restricted possibilities to manipulate and de- 
tect material particles. In particular only position mea- 
surements are easily doable. Existing proposals there- 
fore rely either on reduced fluctuations in the center of 
mass and relative motion 0, in the spirit of Einstein, 
Podolsky, and Rosen (EPR) [J] , or on the violation of a 
Bell inequality Q. Both approaches have drawbacks. 
The former is based on correlations that appear invit- 
ingly easy to explain in terms of a classical (nonquan- 
tum) model, and the latter requires interferometers to 
complete the measurements _6i]. If one restricts oneself 
to elementary position measurements, the states violat- 
ing a Bell inequality maximally seem to be hard if not 
impossible to implement experimentally Q- 

In view of the great success and the compelling power 
of single-particle interference experiments, it is natural 
to ask whether, instead of violating a Bell inequality, it 
is experimentally easier to establish entanglement in the 
motion of material particles by means of similarly impres- 
sive nonlocal matter wave interference. We discuss an ex- 
perimentally amenable class of states which provides such 
nonlocal interference. But conceptually, it is not obvious 
a priori that a nonlocal interference pattern — as intu- 
itively convincing as it may be — can indicate entangle- 
ment, thus strictly excluding the possibility to describe 
the correlations in terms of a separable state. While an 
extensive state tomography could also supply such a rig- 
orous proof, it would be advantageous to possess an en- 
tanglement criterion in terms of observables that can be 



directly read off the interference pattern, merely com- 
plemented by measurements of some likewise accessible 
"conjugate" observables. 

In this Letter we provide such a criterion. To be more 
specific, suppose we hold a two-particle state '^{xi,X2) 
which gives rise to a nonlocal interference pattern when 
subjected to joint position measurements. 
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where the envelope w{xi — xq) localizes particle 1 with 
an uncertainty ax ^ X in the vicinity of xq, and simi- 
larly particle 2 around —xq. Obviously, the interference 
pattern describes correlations in the relative coordinate 
a^rci — xi — X2 of the two particles. But are these corre- 
lations necessarily a signature of entanglement? In the 
case of EPR states (e.g., squeezed Gaussian states), en- 
tanglement can be deduced from the reduced fluctuations 
in both the relative coordinate x-^^i and the total momen- 
tum ptot = Pi +P2, since the canonically conjugate oper- 
ator pairs Xj, pj {[xj, pj] = ih), j = 1,2, set lower limits 
to these fluctuations for separable states @, In the 
situation described by ([1]), in contrast, it is not the rela- 
tive coordinate that is "squeezed," but its value modulo 
A. 

We show how this observation can be employed to de- 
rive an entanglement criterion. The key is to identify 
modular variables [l^l as the appropriate pair of con- 
jugate observables. The criterion is rooted in a state- 
independent additive uncertainty relation (UR) for these 
variables, which remedies the problems arising from the 
operator-valued commutator appearing in the Robertson 
UR. We construct a class of non-Gaussian states, denoted 
modular entangled states, which offer natural and robust 
generation protocols and violate this criterion. The in- 
terference pattern in ([T]) is shown to represent only the 
weakest form of nonlocal correlations exhibited by this 
class. 

Multislit interference. — To discuss the prerequisites of 
particle interference and its relation to the modular vari- 
ables it is instructive to recapitulate the single-particle 
case first. Ideally, the (transverse) state immediately 
after passing an aperture of N slits is described by a 
superposition of TV spatially distinct state components 
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where L denotes the sht separation. The particular shape 
of the single-sht wave function is irrelevant for our 
discussion provided its spatial width ax satisfies ax ^ 
L. This guarantees that the envelope of the resulting 
fringe pattern varies slowly on the scale of a single fringe 
period and thus encloses a large number of interference 
fringes. Note that the state ^ can equally be read as a 
longitudinal superposition of comoving wave packets. 

The subsequent dispersive spreading of the N wave 
packets during the free propagation to the screen results 
in their overlap and interference, yielding in the far-field 
limit the characteristic interference pattern on the screen. 
In terms of the initial state ([5]), this position measure- 
ment at asymptotic times corresponds to a formal mo- 
mentum measurement of p = m{x— (■0ms|x|''/'ms))/^- The 
probability distribution 



\{p\4'Ms)\'' = \{pm'FM[^ 
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exhibits the fringe pattern Fn{x) ^ l + {2/N) Ylf^ii^- 
j) cos{2TTjx). In case of TV = 2 Eq. ^ reduces to the 
sinusoidal fringe pattern of the double slit, whereas for 
N > 2 one obtains the sharpened main maxima and 
suppressed side maxima characteristic for multislit in- 
terference. This reflects a tradeoff between the num- 
ber of superposed wavepackets N and the uncertainty 
of the phase of the interference pattern, in analogy to 
the tradeoff between the variances of a conjugate vari- 
able pair. A similar tradeoff exists between the number 
of fringes M ss apL/h L/ax covered by the envelope 
of the interference pattern and the width-to-spacing ratio 
ax/L « 1/M. 

Modular variables. — These mutual relationships be- 
tween the multislit state ([2]) and the resulting interfer- 
ence pattern ([3]) are captured best by splitting the po- 
sition (momentum) operator into an inte ger component 
Nx (Np) and a modular component x (p) [l^, 
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where x = {x + i?/2)mod^ — £/2 and p — {p + 
h / 2i)mod{h / i) — h/2£. (For convenience, we define the 
modular variables symmetrically with respect to the ori- 
gin.) Recent applications of the modular variables are 
discussed in [ll| - [l3j . 

For the multislit state ([2]) the adequate choice of the 
partition scale is given by ^ = L. The probability dis- 
tribution ^ can then be written as \{p = Nph/L + 
p|V'Ms)P ~ \{p = Nph/L\i(;)\'^FN{pL/h), which indicates 
that the modular variables isolate different characteris- 
tic aspects of interference: the periodic fringe pattern is 
described by the modular momentum p, its envelope by 



the integer momentum Np. Similarly, Nx describes the 
distribution of wave packets in ^ and x their (common) 
shape. 

The modular variables x, p have the remarkable prop- 
erty that they commute, [x, p] = 0, despite originating 
from conjugate observables [10) [3 • The common eigen- 
states \x,p) of x and p with eigenvalues x and p read 
\x,p) = y^X^riez *5xp(ipn£/?i)|n^ -I- or, equiv- 

alently, \x,p) = y/T/iexp{-ipx/h)J2,nez\'"^^/^ + p)p- 
The tradeoff between the number of superposed wave 
packets N and the phase of the interference pattern is 
now reflected by a conjugate relationship between the 
integer position N^; and the modular momentum p. 
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Similarly, the tradeoff between the width-to-spacing ratio 
ax/L and the number of covered fringes is described by 
the commutator of the modular position x and the integer 
momentum 



[x, Np] = 
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The projection operators on the right-hand side of ([5]) 
and © result from the boundedness of the modular vari- 
ables and are indispensible to ensure the validity of the 
Robertson UR. This is similar to the relationship between 
an angular position operator and its conjugate angular 
momentum [15|. 

Squeezed modular position states. — The multislit states 
([2|) display their interference in momentum. In view of 
the symmetry between the two pairs (N^,, p) and (x, Np), 
one can construct another class of states where the mod- 
ular variables exchange their roles. This is achieved by 
superposing wave packets that are distinct in momentum 
(instead of position). 
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where {x\ipxa,po) — (2; — xq) exp[ipo(a; — xo)/h] denotes 
a (well-behaved) wave packet that is localized in phase 
space around {xo,po). Nq represents an arbitrary base 
integer momentum. Distinctness of the wave packets re- 
quires that their momentum width ap is smaller than 
their separation in momentum space, ap <C h/X, or, 
equivalently, ax ^ X. A hypothetical "momentum grid" 
with slit width h/d, d ^ X, could, e.g., prepare a state 
with (j3{x) = sine [irx/d] /Vd. 

A position measurement of the states (0 reveals an 
interference pattern with periodicity A, K^ji/'SMp)!^ = 

{x — Xq) pF/v((a;— Xo)/A). Note that xq determines the 
phase of the fringe pattern. Increasing N results in the 
formation of sharp main maxima, which is reflected by 
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4 


10 


100 


SliN) 
S2{N) 


0.0 
0.0 


0.61 
0.30 


0.71 
0.46 


0.79 
0.55 


0.92 
0.76 


0.99 
0.96 



TABLE I: Evaluation of the squeezing functions Si{N) and 
S2{N) for several superposition ranks A'". Si{N) and S2{N) 
describe the squeezing of the modular position x in the single- 
particle case (|8} and of the modular relative position xi — X2 
in the two-particle entangled case (|ll[l. respectively. 



the decreasing variance of the modular position variable 
(now with i — X), 

((Ax)2)sMP = ^[l-^i(A^)]. (8) 

In this sense one may denote the states ^ as squeezed 
modular position states. The monotonically increasing 
squeezing function Si{N) = -{l2/TT^)J2fji{~^y{N~ 
j)/Np < 1 (for xo — 0) is evaluated in TableUfor repre- 
sentative N . Notably, in the limit N oo the variance 
vanishes, indicating perfect squeezing. 

Correspondingly, the variance of the integer momen- 
tum Np increases with N, ((ANp)^)sMP = (N^ — 
1)/12. For iV = 1, however, ((ANp)2)sMP van- 
ishes (since the IV'xq.po) localized on the scale of 
the integer momentum), while ((Ax)^)smp remains fi- 
nite according to ([8]). Validity of the Robertson UR, 
((ANp)2)sMp((Ax)2)sMP > |([x,Np])sMp|V4, thus re- 
quires that the projector on the right-hand side of (jH]) 
renders the Robertson UR trivial for = 1. Indeed, 
we find |([x, Np])sMp| = (^/2^) [l - (1 + (-l)^+i)/2Af] , 
which vanishes for A^ = 1 . This irrelevance of the Robert- 
son UR in the case A^ = 1 impedes its employment in 
the separability criterion presented below. Note that 
|([x, Np])sMp| converges towards the canonical constant 
value, while the projector term in ([6]) is still relevant for 
A' = 3. Its alternating structure can be traced back to 
either minima or (side) maxima of the fringe pattern co- 
inciding with X = A/ 2. 

As an advantage of the modular position squeezed 
states ([7]) compared to the modular momentum squeezed 
states (HI), they exhibit interference by immediate posi- 
tion measurements, while to determine the integer mo- 
mentum one must only distinguish the macroscopically 
distinct components |i/'a;o.po); which is easy once they are 
sufficiently separated by free propagation. At the same 
time, as superpositions of different velocities, they are 
genuine matter wave states without photonic analogue. 

Modular entangled states. — We are now prepared to 
move on to entangled states of two material particles. 
Ultimately, we are interested in states that reveal their 
entanglement by a nonlocal interference pattern similar 
to ([T]). To this end, we introduce two-particle modular 
position entangled (MPE) states, which are defined by su- 
perposing correlated pairs of (counterpropagating) wave 



packets of different velocities, 
1 

I'I'mPe) = -y= ^ \'lpxo,(No+n)h/\)l\'4'-xo.-iNo+n)h/x)2- 

V 71 — 

(9) 

Only for clarity we assume that the particles are spa- 
tially separated, positioned at ixg. Moreover, it is clear 
that one could equally define modular momentum en- 
tangled (MME) states. Such states could be generated 
(to good approximation) by the sequential coherent dis- 
sociation of a diatomic molecule [16]. For convenience 
we consider a superposition of product states; correlated 
components \'^xo,po\-xo,-po) would not modify our con- 
clusions, since the latter are based on entangled integer 
momenta Npj-, i.e., distinctive "bulk" properties of the 
particles. This is in contrast to EPR states, where the 
relevant correlations reside in the microscopic fluctua- 
tions of the center of mass and relative variables. Per- 
forming position measurements on each side, the non- 
separable structure of ([9]) gives rise to an interference 
pattern in the relative position x^d, |(a;i, a;2|^MPE)P = 
\(l){xi - Xo) \'^\4'{x2 + Xo) \'^FNiixi-X2)/X) (with xo.i 
^0,2), or, equivalently, to a squeezing in the modular rel- 
ative position Xrci = xi — X2. 

Modular entanglement criterion. — The correlations in 
iroi and the total integer momentum A'p tot — -^p,i+-^p,2 
exhibited by the MPE states © can be exploited to 
demonstrate the underlying entanglement. In analogy 
to 3, we consider the sum of variances, ((ANp^tot)^)p + 
((Axi.oi)^)p/^^. Using the Cauchy-Schwarz inequality 
one can show that a separable state of motion, p = 

Y.iPiPu®P2i, implies ((ANp,tot)^)p + ((Axroi)^)p/^^ > 
E,,,p.{((ANp,,)2), + ((Ax,)2),/^2}^ ^ith j = 1,2. In 

contrast to ^1, we cannot use the Robertson UR to es- 
timate the remaining sums of variances, since the ex- 
pectation value of the state-dependent commutator ^ 
vanishes when evaluated for an MPE state (jH]). How- 
ever, one can establish a state-independent additive un- 
certainty relation for the modular variables, ((ANpj)^) + 
((Axj)^)/£^ ^ CNp,x > 0. Using this, we immediately get 
the desired criterion, 

((ANp^tot)')p + ^((Ax,el)2)p ^ 2CN„X, (10) 

which must be satisfied by any separable state. The 
constant cn^^x is given by the smallest eigenvalue //q of 
the operator kj ~ N^ ^ -|-x|/£^. The corresponding dif- 
ferential equation in the common eigenbasis of Xj and 
is solved by ipixj^pj) = exp(— 7ra;^/£^)M(— 7r/x/2 -|- 

l/'^,l/2,2T:xyi'^)x{pj), with M{a,b]x) the Kummer 
function and x{P) arbitrary. Continuity requires a van- 
ishing first derivative at Xj = f/2, which implicitly de- 
termines the (discrete) spectrum {^u} of Aj. Its smallest 
eigenvalue evaluates numerically as cn^.x = 0.078235; sec- 
ond order perturbation theory (with Aj expressed in the 
common eigenbasis of Np^-, p^) yields a reasonable an- 
alytic approximation, cn^.x ~ 1/12(1 — 1/15). We note 
that a criterion similar to can be established for 
^x,tot and p,.ci. 
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The MPE states ^ violate the separabihty criterion 
([TU)) for any N ^ 2. Indeed, the resuhing variances read 
((ANp,tot)^)MPE = and 

((Ax,el)'>MPE- y[l-52(iV)], (11) 

(again with £ = X) where the monotonically increasing 
squeezing function 5*2 (TV) = {6/Tr^)J2fSj^iN-j)/Nf < 
1 is evaluated in Table U for several representative N. 
This proves the possibility to deduce entanglement from 
a nonlocal interference pattern. Again, one can achieve 
perfect squeezing in the limit iV — >■ oo; the interference 
pattern ([T]) corresponds to iV = 2. 

The MPE states ^ (and MME states alike) generalize 
single-particle interferometric schemes such as double-slit 
or grid experiments to the case of two entangled parti- 
cles. Aside from the additional requirement to provide 
the correlations between the particles, the MPE states 
thus inherit both the advantages and the challenges of 
such schemes. Similar to any interference experiment, the 
phase xq of the superposed components \ipxo,po) must be 
well controlled, and also all components should share the 
same shape (j){x) [see ([7|)]. (On the other hand, the partic- 
ular shape is to a large extent irrelevant, which leaves it 
to the experimenter to choose easily producible states.) 
Deviations from these conditions result in a visibility- 
reducing blurring of the fringe pattern and thus in an at- 
tenuation of the squeezing of the modular variable. How- 
ever, a simple robustness check, where the MPE states 
are mixed with merely classically (integer momentum) 
correlated states, reveals that for = 2 a classically 
correlated admixture of up to 79% would sustain the vi- 
olation of the separability criterion, corresponding to a 
fringe visibility of merely 21%. This robustness, which 
even improves with increasing A^, underlines the appro- 
priateness of the separability criterion (jlO[) to capture 
entanglement in spatial interference, and it should leave 
sufhcient freedom to cope with possible experimental lim- 
itations. 



A realistic generation protocol for MPE states would, 
e.g., gradually dissociate an ultracold diatomic Feshbach 
molecule such that subsequent dissociation instants pro- 
duce wave packets with staggered kinetic energies 16| . 
Appropriate dissociation pulses can achieve that all of 
these consecutive wave packets meet simultaneously on 
each side. This constitutes an approximate MPE state, 
where the superposed lipxo.po) then realistically differ by 
different stages of dispersion. We checked for N ~ 2 
and lithium atoms that this dispersion-induced shape dif- 
ference can easily be kept under control with realistic 
parameters, yielding an experimentally resolvable fringe 
pattern with A « 100/im and a visibility of 85%. On 
the other hand, a "grid state preparation" of transversal 
MME states, starting, e.g., with an EPR correlated par- 
ticle pair and then each particle passing a grating, would 
provide the identity in phase and shape of the \ipxo,po) 
for any N by means of the state preparation. 

Conclusion. — We presented a scheme to provide and 
detect entanglement in the motion of two free material 
particles. Elementary position measurements at macro- 
scopically distinct sites give rise to a nonlocal interference 
pattern; the nonseparability then follows from reduced 
fluctuations in adapted modular variables. In this sense, 
the scheme allows one to "deduce entanglement from in- 
terference" , and hence to illustrate the wave-particle du- 
ality on a new level including quantum mechanical non- 
locality. 

We emphasize that the modular variables are merely 
a matter of interpretation in our scheme and can be de- 
duced from ordinary position and momentum measure- 
ments. Finally, it is clear that the entanglement crite- 
rion is applicable to any bipartite continuous variable 
system with conjugate operator pairs, e.g., quadrature 
amplitudes of field modes, and could thus offer a valuable 
alternative to existing entanglement detection schemes. 
Homodyning entangled coherent states [13] niay serve as 
an immediate example. 
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